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The successful background-independent quantization of Loop Quantum Gravity relies on the key 
observation that classical General Relativity can be cast into the connection-dynamical formalism 
with the structure group of SU(2). Due to this particular formalism, Loop Quantum Gravity was 
generally considered as a quantization scheme that applies only to General Relativity. However, we 
will show that the nonperturbative quantization procedure of Loop Quantum Gravity can be ex- 
tended to a rather general class of metric theories of gravity, which have received increased attention 
recently due to motivations coming form cosmology and astrophysics. In particular, we will first in- 
troduce how to reformulate the 4-dimensional metric f(R) theories of gravity, as well as Brans-Dicke 
theory, into connection-dynamical formalism with real 517(2) connections as configuration variables. 
Through these formalisms, we then outline the nonpertubative canonical quantization of the f(R) 
theories and Brans-Dicke theory by extending the loop quantization scheme of General Relativity. 
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I. INTRODUCTION 



q-( In the recent 25 years, loop quantum gravity (LQG), a background independent approach to quantize general 
relativity (GR), has been widely investigated [l|-y . It is remarkable that, as a non-renormalizable theory, GR can be 
^Jy non-perturbatively quantized by the loop quantization procedure. This background-independent quantization relies on 
—i , the key observation that classical GR can be cast into the connection-dynamical formalism with the structure group of 
Thus one is naturally led to ask the following questions. Does the non-perturbative quantization scheme 
of LQG apply only to GR? What is the applicable scope of LQG? On this Loops 11 Conference, we just learned from 
the talks by Thiemann et al that LQG is applicable to GR in arbitrary dimensions, since higher dimensional GR could 
also be cast into connection dynamics with certain compact structure group. We are going to answer the question: 
, whether GR is a unique relativistic theory of gravity with the desired connection-dynamical character? 

In fact, modified gravity theories have recently received increased attention due to motivations coming form cosmol- 
^vq ogy and astrophysics. A series of independent observations, including type la supernova, weak lens, cosmic microwave 
■ background anisotropy, baryon oscillation, etc, implied that our universe is currently undergoing a period of acceler- 
i— I | ated expansion 7]. This result has carried the "dark energy" problem in GR. A large number of phenomenological 
models for dark energy have been proposed, such as quintessence, phantom, Chaplygin gas, K-essence etc. In all these 
models, certain scalar fields are added by hand, and their origins are not understood. Moreover, there are indirect 
evidences suggesting that the bulk of the matter of the universe is invisible or dark (see e.g. Q). The strongest 
evidence comes from the rotational velocity of the isolated stars or hydrogen cloud on the outskirts of galaxies. Al- 
though a few candidates for dark matter are proposed, such as sterile neutrinos, neutralinos etc, none of them are 
satisfying. To explain the accelerated expansion of the universe, as well as dark matter, from fundamental physics is 
now a great challenge. It is reasonable to consider the possibility that GR is not a valid theory of gravity on a galactic 
or cosmological scale. Historically, Einstein's GR is the simplest relativistic theory of gravity with correct Newtonian 
limit. It is worth pursuing all alternatives, which provide a high chance to new physics. 

A large variety of models of f{TZ) modified gravity have been proposed to explain the "dark universe" without 
recourse to dark energy and dark matter [9j. The action of metric f(7Z) theories reads: 

%] = \j d 4 xV^f(n) (i.i) 

where / is a general function of the scalar curvature TZ, and we set 8irG — 1. Besides f(7V) theories, a well-known 
competing relativistic theory of gravity was proposed by Brans and Dicke in 1961 10;], which is apparently compatible 
with Mach's principle. To represent a varying "gravitational constant", a scalar field is non-minimally coupled to the 
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metric as 

%,</)] = J d^V^dm-^id^d^). (i.2) 

The scalar field in Brans-Dicke theory (BDT) of gravity is then expected to account for "dark energy". It can 
naturally lead to cosmological acceleration if certain potential term of the scalar is added in the original action (|1.2[) . 
It is also possible to account for the dark matter problem. Furthermore, a large part of the non-trivial tests on gravity 
theory is related to Einstein's equivalence principle (EEP) [ll|, which contains the following three elements, (i) The 
trajectory of a freely falling "test" body is independent of its internal structure and composition, (ii) The outcome 
of any non-gravitational experiments is independent of the velocity of the freely falling reference frame in which it is 
performed, (iii) The outcome of any non-gravitational experiments is independent of where and when in the universe 
it is performed. There exist many local experiments in solar-system supporting EEP, which implies the metric theories 
of gravity. In a metric theory of gravity, spacetime is endowed with a metric, and the trajectories of freely falling 
"test" bodies are geodesies of that metric. Moreover, in local freely falling reference frames, the non-gravitational laws 
of physics are those written in the language of Special Relativity. Actually, besides GR, both f(TZ) theories and BDT 
belong to metric theories of gravity. For metric theories, gravity is still geometry with diffeomorphism invariance as 
in GR. The differences between them are just reflected in dynamical equations and additional variables. Hence, a 
background-independent and non-perturbative quantization for metric theories of gravity is preferable. Since metric 
f(lZ) theories, as well as BDT, are a class of representative metric theories, which have been received most attention, 
we will take them as examples to carry out the extension of LQG to metric theories. Throughout the paper, we use 
Greek alphabet for spacetime indices, Latin alphabet a,b,c,..., for spatial indices, and i,j,k,..., for internal indices. 



II. CONNECTION DYNAMICS OF f(K) AND BRANS-DICKE THEORIES 

By introducing an independent variable s and a Lagrange multiplier cf>, an action equivalent to (jl.ip of f(JZ) theories 
is proposed as 

S[g, h*] = \j <?xV=g(f(8) - <f>(s - K)). (2.1) 

The variation of (|2.1j) with respect to s yields <f> = df(s)/ds = f'(s). Assuming s could be resolved from the above 
equation, action (|2.1[) is reduced to 

S[9, 4>] = lf d 4 xV^m - £(</>)) = I d A xC (2.2) 

where = <fis — f(s). Comparing action (|2.2[) with (|1.2j) . it is obvious that metric f(TZ) theories could be regarded 
as a particular kind of generalized BDT with certain potentials o f <j> and a vanishing coupling parameter oj. The 
virtue of (|2.2p is that it admits a treatable Hamiltonian analysis Qjj]- By doing 3+1 decomposition and Legendre 
transformation, the Hamiltonian of metric f(lZ) gravity can be derived as a liner combination of first-class constraints 
as: 

Hiatal = f N a V a + NH, 

where N and N a are the lapse function and shift vector respectively, and the diffeomorphism and Hamiltonian 
constraints read 

V a = -2D b (p ab )+nd a ^ 

H = ^={ Pahpab ~ ^ + W ~ -P7T) + ±VKm -4>R + 2D a D^), 
y/h <? 6 3 2 

where p ab and it are the momentum respectively conjugate to the spatial 3-metric h a t, and scalar field 4> on the spatial 
manifold E. Although the above Hamiltonian analysis is started with the action (|2.2|) where a non-minimally coupled 
scalar field is introduced, we can show that the resulted Hamiltonian formalism is equivalent to the Lagrangian 
formalism. The symplectic structure reads 



{hab{x),pr d b)} = 8i°8*>P(x,y), {4>(x)My)} = S 3 (x,y). 



(2.3) 
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To achieve the connection dynamics of metric f(7Z) modified gravity, we let 

where K ab is the extrinsic curvature of E, and introduce 

E? = Vhe?, Ki=K ab e bj , 

where ef is the triad such that h a i,efe b = &y. Now we first extend the phase space of geometry to the space consisting 

of pairs (Ef,Kl). It is then easy to see that the symplectic structure (|2.3p can be derived from the following Poisson 
brackets: 

{E?(x),E b (y)} = {Ki(x),K*(y)} = 0, {Ki(x) , E b (y)} = 5%8{x,y). (2.4) 

Thus there is a symplectic reduction from the extended phase space to the original one. Since K ab = K ba , we have 
an additional constraint: 

G jk ee K a[j E a k] = 0. (2.5) 
For second step, we make a canonical transformation by defining 

K = K + 7K, 

where T'^ is the spin connection determined by Ef and 7 is a nonzero real number. Then the Poisson brackets among 
the new variables read 

{Ai(x),E b k (y)} = i8 b Ji5( X ,y), {^(*),^ b (y)} = 0. 

Now, the phase space consists of conjugate pairs (A l a ,E b ) and ((f), tt). Combining Eq. (|2.5[) with the compatibility 
condition: d a Ef + e ijk T{E ak = 0, we obtain the standard Gaussian constraint 

Qi = V a E? ee d a Ef + e ljk AiE ak , (2.6) 

which justifies A l a as an Ashtekar-Barbero su(2)-connection. Had we let 7 = ±i, the (anti-)self-dual complex connec- 
tion formalism would be obtained. The original diffeomorphism constraint can be expressed in terms of new variables 
up to Gaussian constraint as 

V a = -Fl b E b + ird a <P, 
7 

where F* b = 2d[ a A b ^ + t\iA k a A\ is the curvature of A l a . The original Hamiltonian constraint can be written up to 
Gaussian constraint as 

IT — ^lipj (J2 1 1 \~ T>m T ^n] £ jklE k E b 
tlfR ~ 2^afc~(7 + -^) s jmn^a I ~J= 

Similar to GR, the constraints are of first class. Thus the total Hamiltonian is a linear combination of the constraints. 
We have obtained the real sw(2)-connection dynamical formalism of Lorentz f(7Z) gravity fl3l [l4| . Note that a 
connection dynamical formalism of Euclidean f(7Z) theories in Einstein frame was derived in [15| 

We now turn to the BDT. By doing 3+1 decomposition and Legendre transformation, the Hamiltonian of BDT can 
be derived from the action (| 1 . 2|) as a liner combination of first-class constraints as H to tai = / s N a V a + NH, where the 
diffeomor phis m constraint V a has the same form as that of metric f(1Z) theories, while the Hamiltonian constraint 
reads pj, [13] 
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In comparison with metric f(1Z) theories, there is a new kinetic term ^^(D a (/))D a 4> of the scalar field in the above 
Hamiltonian constraint. However, this kinetic term will not affect the canonical transformations to connection dy- 
namics. Following the same canonical transformations in f(lZ) theories, we obtain new conjugate variables for the 
geometry as: A l a = T l a + ^/K l a and Ef — Vhef, as well as the Gaussian constraint (I2.6[) . In terms of the new variables, 
the Hamiltonian constraint reads 

Hun = \ [fI - ( 7 2 + ±)e 3mn KTKz) 

Thus both f(lZ) theories and BDT have been cast into connection dynamics with su(2)-connection as one of config- 
uration variables. 



III. LOOP QUANTIZATION OF f(TZ) AND BRANS-DICKE THEORIES 

Based on the connection dynamical formalisms, the nonperturbative loop quantization procedure can be straight- 
forwardly extended to both f(7l) theories and BDT. Since the scalar held in f(7Z) theories and BDT still reflects 
gravity, it is natural to employ the polymer-like representations for both the scalar field and the connection. The 
following quantum kinematical structure is valid for both f(7l) theories and BDT. 

For the geometry sector, we have the unique diffcomorphism and internal gauge invariant representation for the 
quantum holonomy-flux algebra 1 . The kinematical Hilbert space of geometry reads 

with the so-called spin-network basis T a (A) = T„ .i.m.n(A). The spatial geometric operators of LQG, such as the area 
[HI, the volume and the length operators (20|, [2l[ are still valid here. Hence, the important physical result that 
both the area and the volume are discrete at quantum Kinematical level remains true for loop quantum f(lZ) gravity 
as well as loop quantum Brans-Dicke gravity. 

For the polymer- like quantization of the scalar field [22], one extends the space U of smooth scalar fields to the 
quantum configuration space U. A simple element U may be thought as a point holonomy: U\ = exp(i\(j)(x)) at 
point idS, where A is a real number. By GNS structure, there is a natural diffcomorphism invariant measure dfj, on 
U. Thus the kinematical Hilbert space of scalar field reads 

H s k Z=L 2 (U,dn), 

with the scalar-network basis 

XjEX 

where X = {xi, . . . , x n } is an arbitrary given set of finite number of points in E. The total kinematical Hilbert space 
for ,f(lZ) and Brans-Dicke gravity reads Whin ■— Hiin ® w ^h an orthonormal basis 

T a , x (A, cj)) = T a (A)®Tx(^. 

The holonomy h e (A) = V exp J e A a , flux E(S,f) := j s € a bcEf /*, point holonomy U\, and smeared momentum 
tt(R) := f R d 3 xir(x) of scalar field become basic operators in T-Lkin- Their actions read respectively 

h e (AMA,4>) = h e (A)*(A,<t>), E(SJ)*(A,<t>) = ih{E(SJ),V(A,<t>)}, 
U x {(t>{x))^{A,(t>) = exp(i\<Hx))*(A,fl, n(R)^(A,ct>)=ih{w(R),^(A,(f>)}. 

As in LQG, it is straightforward to promote the Gaussian constraint G(A) to a well-defined operator. The kernel 
of (?(A) in Hkin is the internal gauge invariant Hilbert space T-Lg, with gauge invariant spin-scalar-network basis: 



1 See e.g. talks by Sahlmann, Giesel, Lewandowski on Loops 11 Conference in Madrid 
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T s<c {A,ct>)=T s 

={a,j,i){A) ® Tx((f>)- Since the diffeomorphisms of £ act covariantly on the cylindrical functions in Hg, 
the so-called group averaging technique can be employed to solve the diffconiorphism constraint |2|-|4(. Thus we can 
also obtain the desired diffeomorphism and gauge invariant Hilbert space Hoiff for f(TZ) and Brans-Dicke gravity. 

The quantum dynamics is a nontrivial issue. To implement the Hamiltonian constraint at quantum level, we need 
to deal with the metric f(lZ) and Brans-Dicke theories separately. The smeared version of the Hamiltonian constraint 
(|2.7p of metric f(Ti) theories can be written in regular order as: Hfji(N) = X)j=i Hi(N). Comparing it with that 
of GR in connection formalism, the new ingredients that we have to deal with are <fi(x), -1 (a;), £(</)) and the four 
terms: H 3 (N), H±{N), H 5 (N), H 7 (N). By introducing certain small constant Ao, an operator corresponding to the 
scalar (j)(x) can be defined as 

fa) = ^t(^a o (0(x)) - tf-AoWaO)). 

ZIAq 

The ambiguity of Ao is the price that we have to pay in order to represent field cj> in the polymer-like representation. 
To further define an operator corresponding to <fi (x), we can use the classical identity 

r X (x) = sgn[4>] (rhgnlflMx^MR)})^ , (3.1) 

for any rational number I £ (0,1). For example, one may choose I = | for positive <f>(x) and replace the Poisson 
bracket by commutator to define 

4>- 1 (x) = (^[^(x),7r(R)}f. 

Then all the functions £(</>) which can be expanded as powers of 4>{x) have been quantized. For other non-trivial 
types of £(4>), we may replace the argument <fi by <j>, provided that no divergence would arise after the replacement. 
In the case where divergence does appear, there remain the possibilities to employ tricks similar to Eq. (|3.1j) to deal 
with it. Hence it is reasonable to believe that most physically interesting functions £(</>) can be quantized. Then it is 
straight-forward to quantize the term Hq (N) in the Hamiltonian constraint as an operator acting on an basis vector 
T a ,x as 

H 6 (N)-T a>x = ^ N(v)i(4>(v))V v -T a , x . 

v£V(a) 

Moreover, by the regularization techniques developed for the Hamiltonian constraint operators of LQG and polymer- 
like scalar field [111, all the terms H 3 (N), H 4 (N), H 5 (N),H 7 (N) can be quantized as operators acting on cylindrical 
functions in T-Lkin in state-dependent ways. The regularization procedure involves re-expressing some variables by 
Thiemann's trick, point-splitting, triangulating £ in adaptation to some graph a U X underling a cylindrical function, 
and replacing connections by holonomies. Replacing the classical variables and Poisson brackets by the corresponding- 
operators and commutators respectively, we can obtain the regulated Hamitonian constraint operator Hj R (N). It 

turns out that the regulators of H^(N) , Hl(N) , H^(N) can be straightforwardly removed. However, the action of 
fflf(N) on a basis vector T a _x is graph changing. It adds a finite number of vertices with representation Ao at 
t(sj(«)) = e for edges ej(t) starting from each high(> 3)-valent vertex of a. The family of operators Hj{N) fails to 
be weakly convergent when e — > 0. However, due to the diffeomorphism covariant properties of the triangulation, 
the limit operator can be well defined via the so-called uniform Rovelli-Smolin topology induced by diffeomorphism- 
invariant states. The operators corresponding to the four new terms in the Hamiltonian constraint act on a basis 
vector as follows. 

H 3 (N)-T atX = -^i— N(v)r\v)[H E ^l),^ v }[H E ^l),^/¥ v }-T a ,x, 



UN)-T a , x = - £ — *(„) £ e{ SL , SMl s N )e{ Sl , Sj , SK )e LMN 

vGV(a)DX 1 ^ ' V ' v(A)=v(A')=v 

x Tr^^^A)^- 1 ^),!])^^^^)^^), (Vv) 3/4 }h SN (A)[K^ A)1 (Vv) 3/4 })e IJK 
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2 ls N(v) 

H 5 (N)-T a , x = E 35 e( ifl )6 E ( v )2 ir ( v ^( v ^( v ^ E e{si,sj,s K )e(s L ,SM,SN) 

v€V(a)nX ~ ^ ' ^ u(A)=ti(A')=» 

x e"*Tr(/* sl(A) [^ (A) ,^ 

x ^(^.(AofeVo- (^) 1/2 ]^m(A')[^ (A0 , (V v ) 1/2 ]h SN(A/) [h^ {Aiy (K) 1/2 ]) • r a , x , 
^r(iV).T Q)X = - £ 372?A 2 o 7 g;4 ( ^ E * E 

»gV(o) ' v y w 6(e)=u «(A)=o 
X e^E/^OM^CAO)) (^A O (0(* S X(AO)) - ^o(^K(AO))) 

x Tr^V^fe^, (V;) 1/2 ]/i SJf (A)[^ (A) , (K) 1/2 ]) • T a , X) 

where e(sj, sj, sk) : = sgn(det(s/SjSi<- )(v)) takes the values (+1, — 1, 0) if the tangents of the three segments sj, sj, sk 
at d (in that sequence) form a matrix of positive, negative or vanishing determinant, and E(v) = ^"3^ with n(v) 

denoting the valence of the vertex v. The whole Hamiltonian constraint operator H(N) — J2 V ^(N) v is internal 
gauge invariant and hence also well defined in the gauge invariant Hilbert space T-Lg- 

Although H(N) can dually act on the diffeomorphism invariant states, there is no guarantee for the resulted states 
to be still diffeomorphism invariant. Hence it is difficult to define a Hamiltonian constraint operator directly on the 
diffeomorphism invariant Hilbert space T-LDiff- The idea of resolution is to introduce the master constraint as in LQG 

M :— - f ^ 



2 7s y/\deth(x)\ 
One then gets the master constraint algebra as a Lie algebra: 

{V(N), V(N')} = V([N,N'}), 
{V(N), M} = 0, 
{M, M} = 0, 

where the subalgebra of diffeomorphism constraints forms an ideal. So it is possible to define a corresponding master 
constraint operator on H.Diff- The master constraint can be regulated via a point-splitting strategy as [25]: 

1 f js.. f * H{x) H(y) 



M e =- d A y d 6 x X e(x - y) 



Introducing a partition V of the 3-manifold E into cells C, we have an operator Hfj a acting on gauge invariant 
spin-scalar-network basis T s<c via a state-dependent triangulation as 
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H Cta -T s , c ={H c>GR + Y J Hh,i)-T s , c = E Xc(v)( E ^t« + E^) ,T «- 

i=3 v£V(a) v{A)=v i=3 



where a denotes the underlying graph of the spin- network state T s . Here, 

32 1 - - x l~~ — 

H E GR,v = ^^ e ( S /'^' S ^) e/JK (^^( 1; ) 71r (V.7(A)^K(A)^; K 1 ( A)' yVu;]) 



3E( 



( ifi )3 7 3 $ 1 ( w ) + "f 2 4>(v))Tr(h Sl (A) [h Si \ A} ,K]h Sj{A) [h s } (A) , K] h SK (A) [fcJ (A) , , 
1 fi 
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and HI v , -fff v , iff v , H^ v can be denned similarly [14| . The family of operators Hq are cylindrically consistent up 

to diffeomorphisms. The inductive limit operator He can be well defined by the uniform Rovelli-Smolin topology. 
Thus we can define master constraint operator M. on diffeomorphism invariant states as 



(M$ Diff )T StC = v ^ e ,^Diff[\ £ H C (H^T S , C }. 



It is obvious that Ai is diffeomorphism invariant. For any given diffeomorphism invariant spin-scalar- network state 
T[ s c ], the norm H-MT^^H^^ is finite. So M. is densely defined in %Diff- Moreover, M. is positive and symmetric 
and hence admits a unique self-adjoint Friedrichs extension. It is then possible to obtain the physical Hilbert space 
of quantum f(lZ) gravity by the direct integral decomposition oil-Luiff with respect to M.. 

Now we turn to the quantum dynamics of BDT. Comparing the Hamiltonian constraint (|1.2[) of Brans-Dicke gravity 
with that of metric f(TZ) gravity in connection formalism, the only new ingredient that we have to deal with is the 
kinetic term: 



H k (N)= [ d 3 xN^-Vh(D a ^)D a q 



By the regularization techniques similar to those for f(lZ) gravity, it can be quantized as a well-defined operator in 

1 ikin 9-S [16;] 

H k -T aiX = E ¥^4WWEW^ ~ 1{V) £ e( Sl ,s JlSK )e( SL , SM , SN ) 

veV(a) ' V V > y ' W v(A)=v(A')=v 
X e iMAr ^ 1 (^ Si (A)))[t>A o (0(t Si (A))) -E/A.(fe l( A)))] 

x Tr( n h SM{A) [hji {AV (V v ) 3/4 }h SN{A) [h^ {AV (K) 3/4 ]) 



x e 



UK 



Since all the other terms in the Hamiltonian constraint are equal to the corresponding terms of metric f(lZ) theories 
up to coefficients, the whole Hamiltonian constraint operator of BDT has been well defined. For the same reason as 
in f(lZ) gravity, we also wish to define a master constraint operator for Brans-Dicke gravity. The "square root" of 
the new term in the regulated master constraint can be quantized as 

2 15 u> 

H k,v = Q4 4/-\ ^2/■fc^4^V x 2 £ ( S -r> S J' S -ff) £ ( s £> s M,SAr) 

v(A)=v{A>)=v 3 4 7 4 (*A ) 2 ^) 4 £(«) 2 
x Tr(nh SMiA) [h; A ) (A) , (V v ) 1/2 }h SN{A) [h^ {A} , (V v ) 3/4 }) 



x e 



UK 



U^(Hs Sl (A>)Wx„m si{ A>))) - Uxo(<t>(Ss l{ A>)))} 

x Tr(nh sAAI) [h^ iAiy (Vv) 1/2 }h SK (A>)[K« { A'y (Vv) 3/i }), 

from which a positive and self-adjoint master constraint operator can be defined on H-Diff ■ Thus the master constraint 
program is also valid for BDT. 



IV. SUMMARY AND OUTLOOK 



Based on the above loop quantization procedure of metric f(lZ) theories and BDT, we can propose a general 
quantization scheme for 4-dimensional metric theories of gravity. The prerequisite is that the theories should have 
well-defined geometrical dynamics, which means a Hamiltonian formalism with 3-metric h a b as one of configuration 
variables and a closed (first-class) constraint algebra (perhaps after solving some second-class constraints). Without 
loss of generality, we assume that the classical phase space of the theory consists of conjugate pairs (h a b,P ab ) and 
((f>A,ft A ), where (f>A could be a rather arbitrary scalar, vector, tensor or spinor field. Then the quantization scheme 
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consists of the following steps, (i) To obtain a connection dynamical formalism, we first enlarge the phase space by 
transforming to the triad formulation as 

(h ab ,p ab ) - (E? = Vhh ab e aj ,Ki = K ab e bj ), 

where K a b — -S^ipab ~ \ph a b) and K a ^E^ — 0. Then we make a canonical transformation to connection formulation 
as: 

(E^Ki)^(E^Ai = Ti+ 7 Ki), 

with the Gaussian constraint, T> a Ef = d a Ef + €ijkA J a E ak = 0, appeared by construction. It is straightforward to 
write all the constraints in terms of the new variables, (ii) For loop quantization, we first find the polymer-like 
representation of the fields (</>a, ft A ), together with the LQG representation of the holonomy-flux algebra. Then the 
kinematical Hilbert space reads Hkin '■= Hkin ® ^-tiw wnere the basic operators and geometrical operators could be 
well defined. By implementing the Gaussian and diffeomorphism constraints as in standard LQG, we could get the 
gauge and diffeomorphism invariant Hilbert space as: Hkin Hg ^ Hmff- To implement quantum dynamics, one 
may first construct the Hamiltonian constraint operator at least in Hg- Usually it could not be well defined on Hmff- 
Then we can construct master constraint operator in HoijS by using the structure of the Hamiltonian operator, (iii) 
One may try to understand the physical Hilbert space by the direct integral decomposition of HDiff with respect to 
the master constraint operator, (iv) One may also do certain semiclassical analysis in order to confirm the classical 
limits of the Hamiltonian and master constraint operators as well as the constraint algebra. The low energy physics 
is also expected in the analysis, (v) Finally, to complement above canonical approach, we can also try the covariant 
path integral (spin foam) quantization. It should be noted that in the present work we only finished steps (i) and (ii) 
for metric f(7V) theories and BDT. 

To summarize, our main results in two folds. First, the 4-dimensional connection dynamics of metric f(TZ) theories, 
as well as Brans-Dicke theory, have been obtained by canonical transformations from their geometrical dynamics. Thus 
GR is not the unique theory of gravity with connection dynamical character. Second, due to the su(2)-connection 
dynamical formalism, the 4-dimensional metric f(7V) theories and Brans-Dicke theory have been nonperturbatively 
quantized by extending LQG scheme. Hence, the non-perturbative loop quantization procedure is not only valid for 
GR but also valid for a general class of metric theories of gravity. In fact, it is not difficult to extend LQG further 
to general scalar-tensor theories of gravity [llj ■ Moreover, since higher dimensional scalar-tensor theories of gravity 
have well-defined Hamiltonian geometrical dynamics, and the symplectic reduction of Bodcndorfer-Thiemann-Thurn 
connection formalism to metric formalism does not depend on dynamics, LQG may also be extended to higher(> 4) 
dimensional scalar-tensor theories of gravity [26| . 

Of course, there are still many open issues on the extension of LQG. It is desirable to find suitable actions for the 
connection dynamics of f(TZ) theories and that of Brans-Dicke theory. We will explore the applications of loop quantum 
f(7Z) and Brans-Dicke theories to cosmology and black holes in future work [27|. It is also desirable to quantize metric 
theories of gravity by the covariant spin foam approach (28|. To conclude, our conservative observation is that LQG 
could be applicable to metric theories of gravity (with well-defined geometrical dynamics) in arbitrary dimensions. 
However, a caution arises from the difference between the 4(or 3)-dimensional and higher-dimensional connection 
dynamical formulations of GR. In particular, in 4-dimensional case, we have both Ashtekar-Barbero connection 
dynamics and Bodendorfer-Thiemann-Thurn connection dynamics of GR. Are the quantum theories corresponding to 
them unitary equivalent to each other 2 ? If the answer is negative, are there any theoretical criteria for judging them? 
Is Bodendorfer-Thiemann-Thurn connection dynamics only preferable in higher dimensions? Another interesting 
question is whether LQG can be extended to non-metric theories, e.g., metric-affine f(lZ) theories. All these open 
issues are fascinating and deserve future investigating. 
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